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ABSTRACT 

We investigate constraints on the distribution of dark matter in the neighbourhood of the Galactic Centre that may eventually be 
attained with the high-resolution Event Horizon Telescope (EHT). The shadow of a black hole in vacuum is used to generate a toy 
model describing how dark matter affects the size of the shadow of the supermassive black hole located at the Galactic Centre. 
Observations by the EHT may constrain the properties of the dark matter distribution in a possible density spike around the black 
hole. Current uncertainties due to both the resolution of the telescope and the analysis of stellar orbits prevent one from discerning 
the effect of dark matter on the measured size of the shadow. The change in the size of the shadow induced by dark matter can be 
seen as an additional uncertainty in any test of general relativity that relies on using the angular size of the shadow to estimate the 
Schwarzschild radius of the black hole. 
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^ 1. Introduction 

o 

• In the framework of the standard model of cosmology, dark mat- 
*^ ter comprises about 23 percent of the energy content of the uni- 
^Hverse. Consequently, unravelling the mystery of the nature of 
O dark matter is one of the most challenging problems of modern 
S— j cosmology. In this context, astronomical observations can pro- 
c/) vide significant constraints on the dark matter distribution. The 
& dawning of a new era of telescopes with high resolution and sen- 
^ J sitivity has widened the scope of precision experiments that may 
be used to set constraints on the properties of dark matter. In 
J> this work, we focus on one such experiment. We consider an 
t-h experiment focusing on the Galactic Centre, as the centre of the 
^sO Milky Way has become a favoured region for searches of dark 
OO matter signatures. The dark matter density is indeed expected 
to be higher in the neighbourhood of the Galactic Centre, ac- 
y—{ cording to N-body simulations. Moreover there is likely to be 
i— H an enhancement of dark matter, if indeed in the form of weakly 
{N| interacting particles, around the supermassive black hole Sgr A* 
t-H located at the centre of the Milky Way. 

The very high and still increasing resolution of the Event 
• i-H Horizon Telescope (EHT) is about to allow one to get a picture 
^ of the shadow of this black hole. Such an imaging survey may 
then be used to constrain the properties of a dark matter spike 
around the centre. This is of especial interest because such a 
geometrical measurement avoids any knowledge of the precise 
nature of the dark matter, such as its decay or annihilation prop- 
erties. One needs a model describing how the dark matter around 
the centre modifies the size of the shadow, in order to be able to 
make quantitative comparisons with future EHT pictures. 

In Sect. 2, we recall the notion of shadow of a black hole in 
vacuum, before presenting in Sect. 3 a toy model that describes 
how the presence of dark matter around a black hole would affect 
the bending of light induced by the compact object, and how 
this could be used to extract constraints on the properties of the 
dark matter distribution. Finally, we analyse the uncertainties 



involved in the process of inferring the properties of the black 
hole from the measurement of the size of the shadow. 



2. Shadow of a Schwarzschild black hole in vacuum 

The shadow of a black hole is closely related to the way that 
light rays are bent by the black hole, and more particularly to the 
deflection angle. Starting from the Schwarzschild metric, which 
describes the geometry of a non-rotating black hole in vacuum, 
one can show the relation between the impact parameter b and 
the closest approach distance r m of a photon (Weinberg 1972): 



(1) 



where rs is the Schwarzschild radius of the black hole. It turns 
out that the deflection angle features a divergence for a particular 
value of the impact parameter. This value defines the shadow of 
the black hole. More precisely, this divergence corresponds to 
light rays infinitely bent by the black hole and thus performing 
an infinite number of loops around it. Such orbits are unstable 
under small perturbations, which results in the photons eventu- 
ally crossing the horizon and falling onto the singularity. There- 
fore, the shadow represents the minimum impact parameter of a 
photon escaping the attraction of the black hole. Conversely it is 
the minimum impact parameter of a photon coming from infin- 
ity for it to end up on the photon sphere (Bozza 2010), as seen 
in Fig. 1. So the shadow turns out to be the effective black disk 
seen by an observer looking at the black hole. The shadow bor- 
der can thus be determined analytically as it corresponds to the 
minimum of the function b(r m ). The photon sphere has a radius 
given by 1.5 r s . This value is obtained by minimizing b, and the 
corresponding value of the impact parameter yields the radius of 
the shadow: r shadow = 3 V3/2 r s « 2.6 r s . 
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Fig. 1. Shadow of a black hole. The radius of the shadow is the 
minimum impact parameter of a light ray escaping the black hole. So 
the shadow is a disk representing the black hole as seen by the observer. 
The circular orbit lies on the so-called photon sphere. The black circle 
represents the horizon. 

Consequently, the shadow is the main observable feature of 
the black hole in a direct imaging survey, and this is what the 
EHT collaboration aims to observe in the near future, using the 
technique of Very Long Baseline Interferometry (VLBI). More 
precisely, the EHT project consists in phasing up millimetre and 
sub-millimetre telescopes scattered over the world and which, 
together, will form an effective high-resolution Earth-sized tele- 
scope (Doeleman 2010). The EHT array currently comprises the 
Submillimeter Telescope (SMT) in Arizona, the Combined Ar- 
ray for Research in Millimeter-wave Astronomy (CARMA) in 
California and three telescopes in Hawaii: the Caltech Submil- 
limeter Observatory (CSO), the James Clerk Maxwell Telescope 
(JCMT) and the Submillimeter Array (SMA) The network 
is soon to be complemented with the Atacama Large Millime- 
ter/submillimeter Array (ALMA) in Chile. The data recorded 
simultaneously by all these telescopes are then processed by a 
dedicated supercomputer. So far, two baselines have been put 
together, namely the CARMA-SMT and the Hawaii-SMT base- 
lines, in order to achieve an angular resolution of 58 microarc- 
seconds (//as) at 230 GHz. However, the angular Schwarzschild 
radius of the supermassive black hole Sgr A* is 10.3 ± 1.24 //as 
and thus the angular diameter of the shadow is of the order of 
53.6 //as. So the current resolution is not yet sufficient to get 
an image of the shadow of the black hole, but it is quite close. 
Therefore, adding more baselines to the EHT array will improve 
the resolution so that the shadow of Sgr A* should be observed 
in the very near future. Moreover, the resolution is better at 
lower wavelengths, therefore future observations are going to be 
carried out at an additional larger frequency, namely 345 GHz. 
Consequently, the addition of ALMA, along with another base- 
line between the Plateau de Bure interferometer in France and 
the South Pole Telescope (SPT), will allow the EHT array to 
achieve a resolution of 15 //as at 345 GHz. 

Imaging the shadow of a black hole will allow one to test 
the prediction of general relativity for the radius of the shadow 
and study accretion flows in the vicinity of the black hole. But 
it will also allow one to set constraints on the properties of a 
hypothetical distribution of dark matter surrounding the super- 

1 see http://www.eventhorizontelescope.org/ 



massive black hole at the Galactic Centre, and this is the object 
of the next section. 



3. Shadow of a black hole surrounded by dark 
matter 

In this section we present a toy model describing how the shadow 
of a black hole is modified if the singularity is embedded in dark 
matter. 

3. 1 . The modified metric 

In order to determine the behaviour of light rays, one must first 
derive the metric describing the geometry of spacetime induced 
by the black hole embedded in matter. The matter distribution 
is assumed to be spherically symmetric for the purpose of sim- 
plification. Let us start from the general form of a static and 
spherically symmetric metric (Weinberg 1972): 

ds 2 = B(r)dt 2 - A(r)dr 2 - r 2 [dG 2 + sin 2 G dip 2 ) (2) 

and let us determine A and B using the Einstein equations. The 
dark matter halo is assumed to be collisionless and at rest. Under 
these assumptions, the first metric coefficient is given by: 

A»-(l-^)" (3) 

r , , 

where m(r) — I A-nr p(r ) dr is the enclosed mass in the 
Jo 

sphere of radius r. The other unknown coefficient is ruled by 
the following differential equation: 

ZT 2m(r) / _ 2m(r)\" 1 

B r 2 \ r ) ' 

3.2. Shadow of a black hole surrounded by a dark matter 
spike 

A dark matter density profile must now be chosen. Here we con- 
sider dark matter distributed according to a density spike around 
the central black hole. More precisely, the galaxy is believed to 
be embedded in a halo with a moderately steep profile (p cc 
where < B < 2) from N-body simulations (Ullio et al. 2001). 
However, some models also predict a dark matter spike at the 
Galactic Centre, that is a steep density profile within the central 
parsec (Gondolo & Silk 1999). The slow growth of the black 
hole at the Galactic Centre, within the halo, may have given rise 
to such a steep density profile for the dark matter component, 
due to accretion of dark matter by the black hole: p cc r~ y where 
2.25 < y < 2.5. If the dark matter halo has a Navarro-Frenk- 
White (NFW) density profile, which goes like p oc r 1 close to 
the centre, then the spike profile turns out to be p cc r~ 1/3 (Gon- 
dolo & Silk 1999). 

In the remainder of this section, we neglect the contribution 
of the halo itself and focus on the gravitational lensing effects of 
the spike, which should be the dominant component of the sys- 
tem, as far as the shadow is concerned. Therefore, the spike is 
assumed to be surrounded by vacuum in the calculations involv- 
ing the shadow. Note that a halo without a spike would have a 
negligible effect on the shadow, due to a less severe enhancement 
of the dark matter density at the centre. 
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Let us consider the following density profile for the system 
comprising the black hole of mass Mbh (equivalent to a point 
mass) and the dark matter spike of mass Mspike and radius R: 



6(r) 

p(r) = M BH ^2 +ar 



(5) 



where y will be taken as 7/3, and a is determined by the normal- 
ization condition that the integral of the spike density be equal to 
the total mass of the spike. Note that the uncertainty on y leads 
to changes in the size of the shadow which are negligible with 
respect to all other uncertainties by several orders of magnitude. 
Therefore with this model the differential equation for B takes 
the form: 



2M R 



B' 
~B 



1 + q 



i - 



2Mb 



1 + q 



(6) 



where the mass ratio is q — M^^/Mbh- Considering that R 
should be much larger than the Schwarzschild radius r$ = 2M BH , 
and q should be smaller than 1, the second term in each bracket 
is going to be a perturbative term for small radii (r ~ r$), rele- 
vant for the derivation of the radius of the shadow. Consequently, 
equation (6) can be approximately integrated using the fact that 
the numerator in the right-hand side is the derivative of the de- 
nominator, except for the small term containing q and R. This 
leads to: 



B*B Q 1 



r 



I r \ 3 ~y 



(7) 



for r < R, where Bq is a constant which can be obtained using 
the expression of the metric coefficient B in vacuum: B(r) = 
1 - rsjat/r for r > R, where rs ;to t = fs(l + q) is the total 
Schwarzschild radius corresponding to M tot = Mbh + Mpike- 
Therefore continuity of the metric at the boundary of the spike 
r — R leads to Bq = 1. So in the following, the metric used to 
describe the geometry induced by a black hole surrounded by a 
dark matter spike will be determined by the coefficients 



A(r) = ' </ 



Bir) « AirY 



(8) 



where the first equation is exact and the second equation holds 
for parameters of the spike such that the perturbative approach 
is valid. More precisely, to make sure that this approximation 
was sensible, we numerically solved equation (6) using a Runge- 
Kutta scheme to compute the actual function Bir), given the 
boundary condition BiR) — 1 - r^, tol IR. Then we computed the 
product of functions Air) given in (8) and Bir) obtained numer- 
ically. It turned out to be equal to 1 with an error of 10~ 3 for 
q = 0.1 and R — 10 3 r$ which are sensible limits from a physical 
point of view (as the spike is predicted to be both large and light 
with respect to the black hole). Such an error will be insignifi- 
cant as the radius of the shadow will be given with a much larger 
uncertainty than one percent (see below). Moreover, these val- 
ues of q and R are the very limits of this perturbative approach 
and for most values of the parameters, the error is actually even 
smaller than 10~ 3 . 

With these approximations, the metric can be derived ana- 
lytically for the density profiles considered here. Therefore, the 



approximations detailed above allow us to pursue the calcula- 
tion and obtain an analytic expression for the size of the shadow 
'"shadow This allows for faster numerical calculations. Then, con- 
straints on the properties of the system can be extracted from 
EHT pictures by measuring the deviation of the actual shadow 
with respect to the expected shadow of a black hole in vacuum. 
Consequently, uncertainties on the measurement of r s h a dow are 
going to play a crucial part. 

First of all, one needs to compute the size of the shadow for 
different values of q and R. Similarly to the case of a black hole 
in vacuum, the impact parameter b as a function of the closest 
approach distance r m is given by: 



b(r m ) 



^/BU~) 



(9) 



Then the radius of the shadow is given by minimizing this func- 
tion with respect to r m . The impact parameter assumes its mini- 
mum for a value of r m such that: 



(mm) - 1 



1 + q 



( (min) \ 



R 



3-n 



(10) 



This equation can be solved perturbatively, considering that 
q (r^^/R) r « 1 for realistic models. We solve it by iter- 
ation, starting from the zeroth order and injecting the solution 
back into the right-hand side of equation (10). Only four itera- 
tions are required to retrieve the solution with a precision higher 
than 10~ 3 . This gives an analytic expression for the radius of the 
shadow r s 



shadow 



b(r ( ™ m \q,R)) as a function of q and R. 
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Fig. 2. Contour plot showing the size of the shadow r shadow as a 
function of the mass ratio q = M sp ^/Mbh and the size of the spike 
R. Solid and dashed lines are the contours defining r shadow to within 
lCT 2 rs, tol . 



3.3. Resulting constraints and uncertainties 

Shown in Fig. 2 are the resulting filled contours of r s h a dow as a 
function of the mass ratio q and the radius of the spike R. It is 
important to note that here r s h a dow is expressed in units of the 
total Schwarzschild radius rs,tot corresponding to the total mass 
M tot = Mbh + M sp -± e = Mbh(1 + <?)■ This is the only quantity 
known from previous studies of the motion of stars at the Galac- 
tic Centre: M tot = 4.3 ± 0.36 x 10 6 M G (Gillessen et al. 2009). 
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This value is in general taken as the black hole mass, but dynam- 
ical studies can not separately measure M sp ± e and Mbh- So we 
consider that the mass given in the literature is the total mass in 
the framework of a black hole surrounded by a dark matter spike. 
Therefore, this allows one to switch from the unknowns R, Mbh 
and M sp ike to R, M tot and q, where M tot is already known. 

Therefore, considering that the size of the shadow obtained 
by solving (10) is expressed in units of the Schwarzschild radius 
r$ of the black hole alone, and as rs tot = r$(l + q), the value of 
/"shadow given by the calculation had to be divided by (1 + q), in 
order to be expressed in terms of rs.tot (which is a known length 
scale, unlike rs). Moreover, expressing Shadow in terms of rs.tot 
allows us to break a degeneracy, as Shadow in units of r$ is very 
close to 2.6 for all values of q and R in the region of interest. 

The aim is then to set constraints on the remaining unknowns 
q and R. However, Fig. 2 shows that there is a degeneracy in 
R. The plot has not been extended to lower values of the ra- 
dius of the spike because R is often predicted to be larger than 
10 3 rs (for instance R ~ 1 pc according to Merritt (2010) and 
rs ~ 10~ 7 pc). Using this degeneracy in R, one can construct a 
more useful graph, displayed in Fig. 3. The quantity of interest 
from the observational point of view turns out to be the angular 
diameter of the shadow Shadow = 2r sha dowARo- 

The black solid line in Fig. 3 is the angular diameter of the 
shadow derived from the previous calculation. This graph pro- 
vides a clearer way of constraining the mass ratio from the value 
of the radius of the shadow. Note that in this study q ranges 
from 10~ 3 to 10 . Going down to smaller values would be 
useless, since r s h a dow(<7) is almost constant below 10~ 3 . More- 
over, the spike is likely to have a perturbative effect with re- 
spect to the black hole, so q should be small enough. That is 
why here q does not go above 10 _1 . Finally, the range consid- 
ered here shows the most compelling variations of the size of 
the shadow: Shadow goes down by 5 //as between 10~ 3 and 10 -1 , 
which sets the minimum uncertainty required to discern an ef- 
fect of dark matter. The yellow shaded area in Fig. 3 denotes 
the uncertainty on the shadow due to the error on the total an- 
gular Schwarzschild radius, Ao^tot = 6.43 //as. The angular size 
of the shadow in this model is indeed derived in terms of the 
total angular Schwarschild radius fls.tot = '"s.tot/^o- This uncer- 
tainty comes from the error on the measurement of the distance 
of Sgr A* from the Sun, R Q = 8.28 ±0.33 kpc, and the total mass, 
M tot = 4.3 ± 0.36 x 10 6 M G (Gillessen etal. 2009). Consequently, 
it turns out that, even before taking account of the resolution of 
the EHT, the uncertainty due to the mass and distance of Sgr A* 
is already larger than the width of the range of values of Shadow 
In order to observe an effect of dark matter, the sum of the uncer- 
tainty on #s,tot an d the error due to resolution should be smaller 
than 5 //as. 

Shown in Fig. 3 are only the best resolutions achievable by 
adding the Hawaii-ALMA and Plateau de Bure-SPT baselines 
to the current baselines (CARMA-SMT and Hawaii-SMT), for 
two observing frequencies at 230 GHz and 345 GHz. Finally, 
the green contour corresponds to a resolution of 1 //as which 
would be ideal to carry out a precise measurement of the size 
of the shadow and observe a possible change induced by dark 
matter. For now a resolution as high as 15 //as, despite being 
huge by astronomical standards, will only allow one to see the 
shape of the shadow. This is the primary goal of the EHT as 
the collaboration aims to shed light on physical processes taking 
place in the vicinity of the black hole, like for instance accretion 
flows. The following step is to measure the size of the shadow 
in order to obtain a direct estimate of the angular Schwarzschild 
radius of the black hole and possibly see an effect of dark matter. 




q=M spike /M BH 

Fig. 3. Angular diameter of the shadow as a function of the mass 
ratio q = M spike /M BH (black solid line). The yellow shaded area is 
the filled contour corresponding to the uncertainty on the total angu- 
lar Schwarzschild radius. The other contours show the uncertainty due 
to resolution, for the baselines (Hawaii-ALMA and Bure-SPT) which 
allow the array to achieve its best resolution, and for two observing fre- 
quencies. The 1 //as resolution is an optimistic case. 

However, with a resolution of 15 //as, the size of the shadow will 
be measured with an error of almost 30 percent. Therefore, with 
such a large error, discerning a deviation from the value of the 
size of the shadow in vacuum is going to prove a challenging 
task. But even measuring the size of the shadow with a high 
enough precision is going to be difficult. As a matter of fact, the 
EHT is going to be able to carry out a direct measurement of the 
total Schwarzschild radius, but with an uncertainty of 15 //as, i.e. 
larger than the error coming from the measurement using stellar 
orbits. So the existing constraints on 0s,tot are not going to be 
improved from the observation of the shadow with a resolution 
of 15 //as. Still, this resolution planned for the near future is 
certainly not the ultimate goal and surely it can be improved even 
more, adding yet longer baselines, with potential EHT sites in 
New Zealand and Africa. 

65 , , — , , — 
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Fig. 4. Angular diameter of the shadow as a function of the mass ratio 
q = M spikc /M BH (black solid line) with only uncertainties coming from 
the error on the total angular Schwarzschild radius (yellow) and an ideal 
resolution of 1 //as (green). 
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Let us now focus on the constraints on the dark matter spike. 
On the one hand, as can be seen in Fig. 4, due to the uncertainty 
on #s,tot an d even f° r a resolution as high as 1 //as, the angular 
size of the shadow can only be known with an error over 7 //as. 
As a result, the error on the size of the shadow is actually too 
large for q to be constrained. The uncertainty on #s,tot turns out 
to be critical and even increasing the resolution of the EHT much 
more would not allow us to set stringent constraints on q. There- 
fore, placing constraints on the dark matter distribution not only 
requires an improved resolution, but also smaller uncertainties 
on the total mass and the distance of Sgr A*. 

On the other hand, due to the degeneracy in R, one can not 
set constraints on the size of the spike directly from this study 
of the shadow. In fact, there is a relation between q and R, but it 
is independent of the calculation of the shadow. So far, for the 
purpose of determining how dark matter affects the shadow, we 
have neglected the smooth distribution (with an NFW profile) 
surrounding the spike, as its effect on the shadow is negligible. 
However, this underlying distribution can be used to derive a 
relation between the mass ratio and the radius of the spike, using 
the continuity of the density at the boundary of the spike: 



\67rp flgal(l + 4)1 

R. I R. \ 2 

with po = p© — - — — + 1 where R gd i = 20 kpc is the ra- 

figal \-Kgal / 

dius characterising the dark matter halo, R G = 8.28 kpc is 
the distance of the solar system from the Galactic Centre and 
Po = 0.3 GeV cirT 3 is the local dark matter density at the Sun. 
Here we have used the fact that R <sc R giii to obtain a simplified 
expression. As q can not be properly constrained at the moment, 
the radius of the spike can not be constrained either, but equation 
(11) provides a way of translating future constraints on q into 
constraints on R, with an uncertainty given by the uncertainties 
on q, M tot and /Jo- 
in any event, although a direct measurement of the effect of 
dark matter is made more difficult by large uncertainties, quan- 
tifying the change in the shadow induced by a dark matter spike 
is useful as it turns out to be an additional source of uncertainty 
when testing the predictions of general relativity by measuring 
the size of the shadow. 

3.4. Comments on some approximations made in this study 

First of all, throughout this work we have neglected the effect of 
stars located around the Galactic Centre on the bending of light 
and more precisely on the size of the shadow of the black hole. 
This seems indeed a sensible approximation since stars around 
the Galactic Centre form a mass distribution that extends much 
further than the dark matter spike. Therefore, even if the lumi- 
nous mass around the Galactic Centre is significant, its effect on 
gravitational lensing is actually negligible due to its being too 
widespread, as can be seen from the previous sections. More- 
over, there is a much less severe enhancement of the density of 
stars at the centre than that expected for dark matter (Merritt 
2010). Stars behave like a background mass distribution, simi- 
larly to the dark matter halo in which the whole galaxy is em- 
bedded, and thus have a negligible contribution as far as changes 
in the size of the shadow are concerned. There is also a central 
nuclear star cluster, which will scatter dark matter particles and 
flatten the density profile to a slope of r~ 3/2 (Gnedin & Primack 
2004). This is important outside a radius where the enclosed 
mass of stars is comparable to the black hole mass, but the dark 



matter profile in the region of interest for the shadow will not be 
affected. 

Then we have made approximations related to the properties 
of the black hole itself. A Schwarzschild black hole is indeed 
an idealized case and from a more realistic point of view, the su- 
permassive black hole at the Galactic Centre should be rotating. 
The effect of the spin of such a Kerr black hole is to compress 
the shadow border on one side and to shift its position (Bozza 
2010). Therefore, it is still possible to compute the radius of the 
shadow, which remains the same as for the Schwarzschild black 
hole for the non-compressed part. Consequently, it was not nec- 
essary to consider the more complicated case of a rotating black 
hole to analyse the effect of a dark matter distribution on the size 
of the shadow. 

Finally, it is worth noticing that this study of the shadow does 
not depend on any assumption about the nature of dark matter. 
The model only depends on a dark matter mass distribution, re- 
gardless of the nature of the dark matter. Therefore, if a positive 
effect on the shadow were to be detected, this may be a sign that 
dark matter is not an artefact due to the way observations are 
interpreted. 

4. Conclusion 

In this work, we have focused on an original way of setting con- 
straints on the properties of dark matter by observing the Galac- 
tic Centre. The study of the shadow of the central black hole 
embedded in a spike of dark matter in principle allows one to 
constrain the properties of the spike, using future observations 
carried out by the Event Horizon Telescope. The EHT is indeed 
close to reaching a high enough resolution to take a picture of the 
shadow. We have presented a model which describes the effect of 
a dark matter spike on the shadow of the black hole Sgr A*. Val- 
ues predicted by this model are then to be compared to the actual 
value of the size of the shadow measured by the EHT, in order to 
discern an effect of dark matter and constrain the mass and the 
radius of the spike. We have quantified the various uncertain- 
ties involved in the problem. It turns out that the uncertainties 
on the mass and distance of Sgr A* which translate into an error 
tot on the shadow, along with the resolution of the EHT array, 
make it impossible to set relevant constraints on the dark matter 
distribution. Consequently, increasing the resolution is not suffi- 
cient to take advantage of a direct observation of the shadow, but 
the focus should also be on reducing the uncertainties on #s,tot- 
This is crucial when trying to discern deviations from the case 
of a black hole in vacuum. Nevertheless, this apparent drawback 
of the method can be turned into an advantage, as this method 
allows one to compute the scale of variation of the size of the 
shadow, 5 //as. This proves to be an additional uncertainty when 
estimating the angular Schwarzschild radius of Sgr A* from the 
size of the shadow, which is one goal of the EHT. 
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